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ABSTRACT 

We investigate the importance of the pseudo-continuum bound-free opacity 
from hydrogen atoms in the atmospheres of cool white dwarfs. This source of 
absorption, when calculated by the occupation probability formalism applied in 
the modeling of white dwarf atmospheres with T e g < 17000 K, dominates all other 
sources of opacity at optical wavelengths. This is unrealistic and not observed. 
On the other hand, a significant flux suppression in the blue part of the spectra 
of cool white dwarfs has been reported, and mainly interpreted as a result of 
the pseudo-continuum absorption from atomic hydrogen. We investigate this 
problem by proposing a new, more realistic approach to calculating this source 
of opacity. We show that this absorption is orders of magnitude smaller than 
that predicted by current methods. Therefore, we rule out the pseudo-continuum 
opacity as a source of the flux deficiency observed in the spectra of cool white 
dwarfs. 

Subject headings: atomic processes - dense matter - stars: atmospheres - stars: 
white dwarfs 



1. Introduction 

Cool white dwarfs are among the oldest objects in the Universe. The detection of large 
populations of these old stars in our Galaxy gives us a unique opportunity to use them 
as chronometers. To achieve this goal, we need reliable atmosphere models for cool white 
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dwarfs. Such models would allow for a better determination of the atmospheric parameters 
of observed white dwarf stars, tracking the evolution of their atmospheric composition, and 
to obtain better boundary conditions for modeling internal structures and the cooling of 
these stars (Fontaine et al. 2000). One problem in modeling cool white dwarf atmospheres 
is the implementation of the Dappen-Anderson-Mihalas (hereafter DAM) model (Dappen et 
al. 1987) in calculating the pseudo-continuum bound-free opacity from the ground state of 
atomic hydrogen. When this source of opacity is calculated following the DAM prescription 
with the occupation probabilities of Hummer & Mihalas (1988, hereafter HM) for white 
dwarfs with T cff < 17000 K, it dominates all other absorption processes at optical wave- 
lengths. Such strong absorption is completely unrealistic and is not observed (Bergeron 
2001; Bergeron et al. 1997). Nonetheless, this process may still be significant in white dwarf 
atmospheres. In fact, Bergeron et al. (1997) suggested that this absorption mechanism may 
be responsible for the observed flux suppression in the blue part of the spectra of cool white 
dwarfs. Following that idea, Bergeron (2001) introduced an arbitrary damping function to 
the DAM pseudo-continuum opacity to obtain good fits of the models to the photometry of 
selected stars. However, the damping function must be adjusted for each star, indicating 
how unsuitable the DAM model is. As a result, the pseudo-continuum bound- free opacity of 
hydrogen atoms is ignored in the modeling of cool white dwarf atmospheres. In this paper 
we readdress this issue by computing this absorption mechanism with a much more realistic 
physical model to assess its importance in cool white dwarf atmospheres. 

The pseudo-continuum bound-free opacity results from the perturbation of hydrogen 
atoms in their ground state by their interaction with other particles in a dense medium 
(Dappen et al. 1987). Such perturbations result in a lowering of the ionization barrier of some 
of the hydrogen atoms and the possibility of a bound-free transition caused by photon with 
energies that are smaller than the ionization potential of the isolated hydrogen atom. DAM 
consider this process in the framework of the occupation probability formalism proposed by 
HM. They constructed a simple model for the optical properties of an interacting hydrogen 
medium. In the weakly ionized atmospheres of cool white dwarfs, the hydrogen atoms 
in higher excited states are perturbed mostly by neutral particles (Bergeron et al. 1991) 
and eventually destroyed (pressure ionized) by excluded volume interactions. The effect is 
stronger for excited states because the size of a hydrogen atom is a monotonically increasing 
function of the principal quantum number. This model for the pseudo-continuum opacity is 
not very precise, as it treats the inter-particle interactions and the associated lowering of the 
ionization potential only in terms of the average sizes of the particles. We will show that a 
quantum mechanical description of the interactions between hydrogen atoms and perturbers, 
absent in the simple DAM model, is necessary to obtain a physically realistic value of the 
pseudo-continuum opacity of hydrogen atoms in the atmospheres of white dwarf stars. 
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Fig. 1 . The most relevant sources of opacity at the photosphere of a pure hydrogen model atmosphere with T e ff = 4600K 

and log g = 7.75 (cgs) - the atmospheric parameters of the white dwarf star LP 380-5 presented on figure 4 of Bergeron (2001). 
The lines represent the pseudo-continuum opacity presented here (red solid), the pseudo-continuum opacity of DAM (blue 
solid), the H — bound-free opacity (dashed), the H2 CIA opacity (long dashed) and the pseudo-continuum opacity derived by 
Bergeron (2001) to fit the spectrum of this star (green solid). 

To accurately account for the effect of lowering the hydrogen ionization barrier in a dense 
medium, we consider short-range, binary collisions. The ionization energy of the perturbed 
hydrogen atom is given by the ionization energy of the H-perturber pair. This energy is 
given by the difference between the exact interaction energy curves calculated quantum 
mechanically of the neutral and singly ionized dimers (interacting pairs). The probability of 
having a pair colliding with a given collision distance can be determined from the interaction 
potential of the neutral dimer. The resulting probabilities of lowered hydrogen ionization 
potential for the physical conditions of white dwarfs atmospheres are orders of magnitude 
smaller than those obtained from the DAM model with the HM occupation probabilities. 

In the next section we discuss the DAM approach to the pseudo-continuum opacity. Our 
new calculations of this absorption mechanism and its application to the physical conditions 
found in the atmospheres of cool white dwarfs are presented in sections 3 and 4, respectively. 
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2. The pseudo-continuum opacity model of Dappen-Anderson-Mihalas 

The DAM model is based on the occupation probability formalism introduced by HM 
for calculating thermodynamical properties of a non-ideal gas. HM introduced the non-ideal 
effects on the equation of state by modifying the internal partition functions of bound species. 
The modified partition function of the hydrogen atom is given by 

Zu = J29me- ei/kBT = J>i7i, (1) 

i i 

where ks is the Boltzmann constant, % indicates the atomic level, gi is the statistical weight, 
€{ is the energy of hydrogen atom counted from the energy of the ground state (i — 0), T 
is the temperature, and uii is the occupation probability. In the HM model, UJi represents 
the strength of non-ideal gas perturbations on an atomic level % . Its value for the phys- 
ical conditions found inside weakly ionized cool white dwarf atmospheres is given by the 
excluded volume interactions model (Bergeron et al. 1991) and decreases for higher principal 
quantum number (HM, Equation 3.4). The number of hydrogen atoms in the excited state 
i is expressed as 

^L = H3L t (2) 

fHot 

where n to t is the total number density of atomic hydrogen. 

DAM extended this model for the calculation of the optical properties of hydrogen 
atoms in the partially ionized plasma. They interpreted the factor 1 — ui{ as the fraction of 
hydrogen atoms with atomic level i being "dissolved," i.e. sufficiently perturbed to describe 
an unbound electron and ion. They assigned this value to the fraction of free-electrons 
states available for a bound-free transition to level i. DAM used this interpretation to derive 
that for the absorption of a photon associated with the transition from level % to level j, 
the probability that level j is bound, and the transition is "bound-bound," is ujj/uji < 1. 
Therefore, the probability of having a "bound-free" transition is 

P bf = 1 - ^. (3) 

To describe the pseudo-continuum absorption to the continuum energy levels localized be- 
tween two discrete but "dissolved" levels of the hydrogen atom, DAM introduced the pseudo- 
level n* with occupation probability uj n *. The latter is calculated by the interpolation be- 
tween the Wj's. Therefore, the probability that, during the % to n* transition, the absorbing 
electron goes into the continuum is given by (Eq. 32 of DAM) 



P b f = l~ — 

UJi 



(4) 
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Fig. A. The interaction energy curves as a function of the inter-particle collision distance for neutral dimers (lower curves) 

and the corresponding single-ionized dimers (upper curves). The different curves for the H — H interaction are for the bound 
(solid line) and anti-bound (dotted line) interaction potentials for H — H and Hj. The different energy curves for the H — H2 
interaction represent the potentials for different orientation of the molecule, with the angle defined between the line connecting 
the pcrturbcr to the center of H2 and the molecular axis. Interaction for the following angles are shown: 90° (solid), 45° 
(dotted), and 0° (dashed). 

The dominant perturbers in the atmospheres of cool white dwarfs are neutral parti- 
cles. Therefore in calculating the DAM pseudo-continuum opacity we used the hard sphere 
model for uji of HM (Section Ilia). In Figure 1 we present the main opacity sources at the 
photosphere of white dwarf star LP 380-5 considered for the potential presence of pseudo- 
continuum opacity by Bergeron (2001). The DAM pseudo-continuum opacity is much too 
high and dominates all other opacity sources by orders of magnitude. This is unrealistic 
and simply not observed. Therefore, the DAM model highly overestimates the probabili- 
ties of pseudo-continuum bound-free absorption (Eq. 4) in the atmospheres of cool white 
dwarf stars. In view of the difficulty with this application of the DAM model of the pseudo- 
continuum opacity, we propose a more realistic approach to calculate this source of absorp- 
tion. 



3. A new model for the pseudo-continuum opacity of H 

The bound-free absorption process results in the ionization of a hydrogen atom. The 
lowering of the hydrogen ionization potential arises from the interaction of the hydrogen 
atom with neighboring particles (Dappen et al. 1987). However, we are interested in the 
pseudo-continuum opacity far from the Lyman edge, i.e. A ^ 1500 A. A significant lowering 
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Fig. 3. The ionization energy of dimers as a function of the inter-particle distance. See Fig. 2 for legend. 

of the hydrogen ionization barrier, by more than ~ 5 eV, is required for bound-free absorption 
at these wavelengths. Such a situation occurs in the case of rare, close collisions for which 
the inter-particle collision distances, r c , are small enough that the probability (P c ) of finding 
a colliding pair with an inter-particle separation smaller than r c is much smaller than unity. 
In this case, multi-particle collisions are insignificant, as the probability of having a close 
iV-particle collision is roughly ~ P^~ l - Therefore, it is sufficient to consider the interaction 
between a H atom and its closest neighbor only. Allard et al. (2004) used this approximation 
to successfully explain the complex shape of the Lyman a line wings detected in the spectra 
of white dwarfs with T er r ~ 12000 K. 

For a given colliding pair the change in the ionization energy results from the formation 
of a temporary dimer, whose ionization energy differs from that of the isolated hydrogen 
atom, Iq = 13.598 eV. This modified ionization barrier I p is given by the ionization energy 
of a dimer calculated at a fixed inter-particle separation r c : 

I P (r c ) = V H +^ pcrt (r c ) - VH_pe r t(r c ), (5) 

where Va- pert (r c ) and l / H+_ pC rt( r c) are the energies of the neutral and singly ionized dimers, 
respectively. This picture is in the spirit of the Franck-Condon principle, which states that 
as a result of large differences between the mass of the absorbing electron and the nuclear 
mass, the nuclei remain fixed during the absorption/emission of a photon (Davydov 1965). 
The differential probability of finding such a dimer with an inter-particle separation between 
r c and r c + dr c in the low density medium (p ^ 0.1g/cm 3 ) is given by (Martynov 1992). 
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dP(r c ) — n pert 




e -V H - peit (r c ,0,4>)/k B T gin r 2 dr< 



(6) 



where n per t is a number density of perturbers, Vk-pert is the interaction energy between 
a hydrogen atom and the perturber localized at the position (r c , 9, (f>) in relation to the 
hydrogen atom, which is assumed to be at the origin of the coordinate system. The dominant 
species, and therefore the main perturbers, in the atmospheres of cool white dwarfs are 
H, H 2 , and He. The sources for the potentials for the interaction of hydrogen atom with 
perturbers are chosen to be: H — H [Kolos & Wolniewicz (1965)], H — H 2 [Broothroyd et al. 
(1991)], and H — He [Shalabi et al. (1998)]. The corresponding potential curves for the singly 
ionized dimers are chosen to be: H^ [Bates & Reid (1968)], H-j~ [Prosmiti et al. (1997)], 
and HeH + [Green et al. (1974)]. These potentials are plotted on Fig. 2. As a result of the 
gerade/ungerade symmetry in the H — H interaction, we have to consider both the bound 
and the anti-bound H — H potential energy curves. For the singly ionized dimers we choose 
the ground state potential energy curves because the upper energy curves (like anti-bound 
state for H^ , Fig. 2.) would result in a much smaller decrease in the ionization energy for 
a given r c . This decrease occurs with much smaller probability than the same change in 
the ionization potential that results from the ionic ground states at a much larger r c , and 
therefore can be neglected. The resulting ionization energy I p (r c ) for a given dimer as a 
function of the collision distance are plotted on Fig. 3. The ionization energy can decrease 
to ~ 5 eV, which occurs at very small r c . Due to the strong repulsion in the neutral dimer 
at short separation (Fig. 2), such a large lowering of the ionization energy occurs with 
very small probability, as for small r c the Boltzmann factor e~ VH - pert ^ kBT << i i n Eq. (6). 
Moreover, in the case of the H — H interaction, only the anti-bound potential leads to a 
significant decrease of the ionization barrier (Fig. 2). 

The probability of a hydrogen atom having an ionization energy I p sufficiently smaller 
for photo-ionization to be caused by a photon of frequency v is 



where the integration is performed over the range of separations r c such that I p {r c ) < hv. 
The resulting probabilities, as a function of photon wavelength A, of having a bound- free 
transition for H for two temperatures and three different compositions are shown in Fig. 
4. The corresponding probabilities from the DAM model (also shown in Fig. 4) are a few 
orders of magnitude larger. We find that for a given density, collisions with H 2 are most 
effective at lowering the ionization energy of H. This is because for a given value of I p (r c ), 
the interaction potential Vh-h 2 ( r c) computed for the orientation angle ~ 90° (see caption of 
Fig 2. for definition) is much less repulsive than Vh-h^c) and Vn-Hc( r c)- 




(7) 
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Fig. 4 . The probability of lowering the bound-free continuum from our model (thick red lines), and the DAM approach 

(thin blue lines), as a function of wavelength for selected temperatures. Calculations for the compositions are shown: (1) 
n H = 10 21 cm -3 (solid line), (2) n H = n. Ha = 10 21 cm -3 (dotted line), and (3) n Ho = 10 23 cm -3 (dashed line). 

The pseudo-continuum opacity is then obtained with 

K bf {v) = T ^-P{v)al f {v) (hu < 13.598 eV) (8) 

where erj^(z/) is the bound-free cross-section of the isolated hydrogen atom extrapolated 
beyond the Lyman edge (DAM), % is the number density of hydrogen atoms, and p is 
the mass density. A proper calculation would consider the bound-free cross-section of the 
dimer as a function of inter-particle separation. To our knowledge this information is not 
available. However, as the photo-ionization cross-sections for H (Palenius et al. 1975) and 
H2 (Ford et al. 1974) differ by less than a factor of 3 and the cross-section is the effective size 
of the absorbing system as seen by the photon, we estimate that the extrapolation of the 
bound-free cross-section of the isolated hydrogen atom beyond Lyman edge introduces an 
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O . The most relevant sources of opacity for representative physical conditions at the photosphere of cool white dwarf 

stars. Upper panels: T = 4000 K, lower panels: T = 5000 K. The lines represent the pseudo-continuum opacity calculated here 
(red solid), the pseudo-continuum opacity of DAM (blue solid), the H — bound-free opacity (dashed), the H2 CIA opacity (long 
dashed), and the He - free- free opacity (dash-dotted). 

uncertainty no larger than an order of magnitude on calculated by Eq. 8. In helium- 

rich atmospheres, where the density can be as high as 2 — 3g/cm 3 , equation (6) should 
be corrected when p ^ 0.1g/cm 3 , by a factor e WH ^ pcrt (Martynov 1992), where wn-pert is a 
thermal potential that arises from the correlations in the dense fluid. Solving the Ornstein- 
Zernike equation in the Percus-Yevick approximation (Martynov 1992) for fluid He, we have 
verified that in the helium-rich atmospheres of cool white dwarfs, e WH ~ pcrt < 100. Therefore, 
due to correlations in the most extreme case, Eq. 8 leads to an underestimate of K^f of 
no more than two orders of magnitude. We will see below that this does not affect our 
conclusions. 



4. Importance of the pseudo-continuum opacity in cool white dwarf 

atmospheres 

Our goal is to investigate the suggestion of Bergeron et al. (1997) and Bergeron (2001) 
that the pseudo-continuum bound-free absorption by atomic hydrogen beyond the Lyman 
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edge may be the missing source of opacity in the atmosphere models of cool white dwarfs 
(T e g ^ 6000 K). For this purpose, we have computed opacities for temperatures, densities 
and He/H composition that are representative of the photospheres of these stars. The 
pseudo-continuum opacity calculated with our model dominates all other sources of opacity 
at wavelengths shorter than ~ 2000 A (Fig. 5). For these cool stars the flux at these 
wavelengths is extremely small and completely negligible (see Fig. 5 of Bergeron (2001)). 
We also note that the uncertainty in the photo-ionization cross-section and our neglect of 
the correlations in dense helium-rich atmospheres do not alter the pseudo-continuum opacity 
enough to make it important at A > 2000 A. 

Fig. 1. reproduces the opacity plot of Bergeron (2001) (Fig. 4 in that paper). This figure 
represents the main sources of opacity for the physical conditions found at the photosphere of 
the white dwarf star LP 380-5. Our pseudo-continuum opacity is several orders of magnitude 
smaller than that necessary to fit the blue spectrum (3000 — 4000 A) of this star (Bergeron 
2001), and is completely insignificant at A > 2000 A. 

5. Conclusions 

The possible presence of an unknown absorption mechanism in the atmospheres of cool 
white dwarfs has been reported by Bergeron (2001) and Bergeron et al. (1997). This opac- 
ity source has been attributed to the pseudo-continuum bound-free opacity from hydrogen 
atoms in their ground state that results from a lowering of the ionization potential because 
of inter-particle interactions in the gas. Opacity models based on the occupation probabil- 
ity formalism highly overestimate this source of absorption in the atmospheres of cool white 
dwarfs. For this reason, it is usually omitted in models. In this paper, we presented a realistic 
model for this absorption mechanism based on binary collisions, exact pair interaction po- 
tentials, and the ionization energy of the colliding pair. We find that the pseudo-continuum 
bound-free opacity decreases very rapidly with wavelength beyond the Lyman edge and be- 
comes completely negligible beyond A ~ 2000 A over the entire range of temperature, density 
and H/He composition relevant to cool white dwarf atmospheres. Therefore, another ab- 
sorption mechanism must be invoked to explain the flux excess of the models in the blue 
part of the spectrum of these stars. 

I thank D. Saumon for useful discussions and comments on this manuscript. This 
research was supported by the United States Department of Energy under contract W-7405- 
ENG-36. 
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